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Protons and helium nuclei are the most abundant components of the cosmic 
radiation. Precise measurements of their fluxes are needed to understand the 
acceleration and subsequent propagation of cosmic rays in the Galaxy. We 
report precision measurements of the proton and helium spectra in the rigid- 
ity range 1 GV-1.2 TV performed by the satellite-borne experiment PAMELA. 
We find that the spectral shapes of these two species are different and can- 
not be well described by a single power law. These data challenge the current 
paradigm of cosmic-ray acceleration in supernova remnants followed by dif- 
fusive propagation in the Galaxy. More complex processes of acceleration and 
propagation of cosmic rays are required to explain the spectral structures ob- 
served in our data. 

Since the discovery of cosmic rays, various mechanisms have been proposed to explain the 
acceleration of particles to relativistic energies and their subsequent propagation in the Galaxy. 
It was pointed out long ago [e.g. (1, 2)] that supemovae fulfill the power requirement to energize 
galactic cosmic rays. Subsequently, models were put forward explaining the acceleration of 
cosmic ray particles via diffusive shock acceleration produced by SN (supernova) shock waves 
propagating in the interstellar medium [see (3) for a review]. 

At the end of the acceleration phase, particles are injected into the interstellar medium where 
they propagate, diffusing through the turbulent galactic magnetic fields. Nowadays, this propa- 
gation is well described by solving numerically (4) or analytically (5, 6) the transport equations 
for particle diffusion in the Galaxy. The galactic magnetic fields mask the arrival direction of 
charged particles, making the cosmic-ray flux isotropic although there are hints of anisotropy in 
the 10-100 TeV range (7). 
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Recent PAMELA measurements of the antiparticle component of the cosmic radiation (8- 
10) have prompted a re-evaluation of possible contributions from additional galactic sources, 
either of astrophysical [e.g. pulsars {11)] or exotic [e.g. dark matter {12, 13)] origin. Detailed 
knowledge of cosmic ray spectra is needed to: a) identify sources and acceleration/propagation 
mechanisms of cosmic rays; b) estimate the production of secondary particles, such as positrons 
and antiprotons, in order to disentangle the secondary particle component from possible exotic 
sources; c) estimate the particle flux in the geomagnetic field and in Earth's atmosphere for 
in-orbit dose estimations and to derive the atmospheric muon and neutrino flux, respectively. 

We present absolute cosmic ray proton and helium spectra in the rigidity interval between 
1 GV and 1.2 TV (Fig. 1, Tables SI, S2), based on data gathered between 2006-2008 with 
PAMELA, a detector orbiting the Earth in a 350-610 km, 70° inclination orbit as part of the 
Russian Resurs-DKl spacecraft {23). 

Our results are consistent with those of other experiments (Fig. 1), considering the statistical 
and systematic uncertainties of the various experiments. There are differences at low (< 30 
GeV) energies caused by solar modulation effects (PAMELA was operating during a period 
of minimum solar activity with solar modulation parameter, 0, of 450-550 MV in the spherical 
force field approximation {24)). PAMELA results overlap with ATIC-2 data {18) between ~ 200 
GV and ~ 1200 GV, but differ both in shape and absolute normalisation at lower energies. The 
extrapolation to higher energy of the PAMELA fluxes suggest a broad agreement with those 
published by CREAM {20) and JACEE (22) but are higher than the RUNJOB {21) helium data. 

To gain a better understanding of the spectra, we have analysed our results in terms of rigid- 
ity instead of kinetic energy per nucleon (Fig. 2 and Tables S3, S4). Two important conclusions 
can be drawn from the PAMELA data. 

Firstly, the proton and helium spectra {J{R)) have different spectral shapes. If a single 
power law, J{R) = AR'^''^ , is fit to the data between 30 GV (above the influence of solar 
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Figure 1 : Proton and helium absolute fluxes measured by PAMELA above 1 GeV/n, compared 
with a few of the previous measurements (14-22). All previous measurements but one (17) 
come from balloon-borne experiments. Previous data up to few hundred GeV/n were collected 
by magnetic spectrometer experiments (14-17, 19) while higher energy data come from calori- 
metric measurements. PAMELA data cover the energy range 1 GeV -1.2 TeV (1-600 GeV/n 
for He). The fluxes are expressed in terms of kinetic energy per nucleon, converted from the 
rigidity measured in the tracker and neglecting any contribution from less abundant deuterium 
(d/p ~ 1%) and ^He (^He/^He ~ 10%). Pure proton and "^He samples are therefore assumed. 
Error bars are statistical, the shaded area represents the estimated systematic uncertainty. 

modulation) and 1 .2 TV, the resulting spectral indices are: 

73^1000 GV,p = 2.820 ± 0.003{stat) ± 0m5{syst), 
73^1000 GV,Hc = 2.732 ± 0.005(stat)t°:[J°^(syst), 
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Figure 2: Proton (top points) and helium (bottom points) data measured by PAMELA in the 
rigidity range 1 GV - 1.2 TV. The shaded area represents the estimated systematic uncertainty. 
The lines represent the fit with a single power law and the Galprop (25) and Zatsepin (2(5) 
models. Details of the models are presented in Tables SI, S2. 



which establishes that there is a significant difference between the two spectral indices in this 
rigidity region. These effects are also seen in Fig. 3 (and in Table S5), where the proton-to- 
helium flux ratio is shown as a function of rigidity. Presenting the results as a ratio reduces the 
possible impact of systematic errors because a number of instrumental effects cancel in the ratio, 
e.g. the estimation of live time and the error associated with the alignment of the tracker and the 
track reconstruction algorithm. The proton-to-helium flux ratio shows a continuous and smooth 
decrease as the rigidity increases. The same ratio cast in terms of kinetic energy per nucleon 
or total kinetic energy exhibits more irregular behaviour (Fig. SI). By applying a power law 
approximation to the two spectra, the ratio can be used to determine the difference between the 
two spectral indices with a smaller associated systematic error, A^r = 7^ — 7|^e = 0.101 ± 
0.0014(stat) ± 0.0001{sys). The ratio is well described by a power law down to rigidities as 
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Figure 3: Ratio of the flux between proton and helium data of PAMELA vs. Rigidity. The 
shaded area represents the estimated systematic uncertainty. Lines show the fit using one single 
power law (describing the difference of the two spectral indices), the Galprop (25) and Zatsepin 
models with the original values of the paper (26) and fitted to the data. Details of the models 
are presented in Tables SI, S2. 

low as 5 GV (green line in Fig. 3). For rigidities R >> (j), the ratio of the two species is 
independent of the solar modulation parameter and allows for the interstellar spectrum to 
be measured in the rigidity range 5-30 GV, where solar modulation effects dominate. Previous 
measurements (14-17, 19) did not have the statistical and systematic precision to demonstrate 
this decrease in the ratio. 

Secondly, as seen in Fig. 4, the PAMELA data show clear deviations from a single power 
law model: 

The spectrum of protons gradually softens in the rigidity range 30-230 GV. In the rigid- 
ity range 30-80 GV, 7|)_8ogv,p = 2-801 ± 0.007 {stat) ± 0.002 {syst), which is lower than 
the value fitted between 80-230 GV: 7|^o-23ogv,p = 2.850 ± 0.015(stat) ± 0.004(s?/st). In 
the case of helium, tIo-sogvhc = 2.71 ± 0.01 {stat) ± 0.002 (syst), which is lower than 
780-230 GV,He = 2-77 ± 0.03{stat) ± 0.004:{syst). We applied Fisher's and Student's t-tests to 
the single power law hypothesis in the range 30-230 GV for both protons and helium (see Sec- 
tion 5 of the Supporting Online Material (SOM (27)) for details). This hypothesis is rejected at 





Figure 4: Proton (left panel) and helium (right panel) spectra in the range 10 GV - 1.2 TV. 
The grey shaded area represents the estimated systematic uncertainty, the pink shaded area 
represents the contribution due to tracker alignment. The straight (green) lines represent fits 
with a single power law in the rigidity range 30 GV - 240 GV. The red curves represent the fit 
with a rigidity dependent power law (30-240 GV) and with a single power law above 240 GV. 



95% confidence level (C.L.). Considering the same rigidity interval in terms of kinetic energy 
per nucleon the Fisher's and Student's t-tests reject a single power law hypothesis at 99.7% 
C.L.. 

At 230-240 GV the proton and helium data exhibit an abrupt spectral hardening. Applying 
Fisher's test and Student's t-test to the proton spectrum above 80 GV, the single power law 
hypothesis is rejected at 99.7% C.L. if only statistical errors are considered. A similar result is 
obtained if the fluxes are increased in line with the systematic uncertainties. If the fluxes are 
instead decreased, the single power law hypothesis is rejected at 95% C.L.. The hardening of 
the proton spectrum occurs at 232I3Q GV with change of spectral index from 78o-232GV,p — 
2.85 ± 0.015 ± 0.004 to 7^232GV,p = 2.67 ± 0.03 ± 0.05. For the helium data, the single 
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power law hypothesis is rejected at 95% C.L. with spectral hardening setting in at 243l3{ 
and a corresponding change of spectral index of 78o-24ogv He = 

2.766 ± 0.01 ± 0.027 and 

7>243GV He = 2.477 ± 0.06 ± 0.03. As a consistency check, we repeated this analysis with the 
three highest energy data points excluded: no changes in the proton and helium results were 
observed. We obtained similar results when using alternative statistical methods such as the 
cumulative sum test (see Section 5.4 in the SOM (27)). 

One of the most striking features of the cosmic rays prior to PAMELA observations was 
their apparently featureless energy spectra. Until now, single power laws, as predicted by the 
shock diffusion acceleration model and diffusive propagation in the Galaxy [see (28) for a recent 
review], could reproduce spectra using similar spectral indices (a fit to the experimental data 
yields 7 ~ 2.7) for protons and heavier nuclei up to energies of about ^ 10^^ eV (the so-called 
'knee' region). Such assumptions are routinely incorporated into common used propagation 
models, such as GALPROP (4), which is widely considered to be the standard model of cosmic- 
ray acceleration and propagation. Our results challenge this scenario (29). As it can be seen 
in Figs. 2 and 3 the GALPROP calculation does not reproduce PAMELA data across the full 
rigidity region. Moreover it is difficult, even with recent models of non-linear shock acceleration 
[e.g. (30, 31)], to produce significant differences in the proton and helium spectra as low as a 
few tens of GV. 

The hardening in the spectra observed by PAMELA around 200 GV could be interpreted 
as an indication of different populations of cosmic ray sources. As an example of a multi- 
source model. Fig. 2 shows a comparison with a calculation (blue curves) by Zatsepin and 
Sokolskaya (26), which was put forward to explain ATIC-2 data (18) and considered novae 
stars and explosions in superbubbles as additional cosmic-ray sources. The parameters of the 
model were fitted to match ATIC-2 data and, consequently, are in disagreement with PAMELA 
data in absolute fluxes and the ratio. If the parameters of this model are fitted to the PAMELA 
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data the agreement can be greatly improved (red curves of Fig. 2 and 3). CREAM also reported 
a direct measurement, albeit with a low statistical and systematic significance, of a change of 
the slope for nuclei (Z > 3) at 200 GeV/n, i.e. at a higher rigidity 400 GV) than our observed 
break in helium spectrum. 

An indication that proton and helium have different spectral indices at high energy (~ 10 
TeV) was reported by JACEE (22). More recently CREAM (20) indirectly inferred [using also 
AMS (77) and BESS (32) data)] that spectral deformation should occur at about 200 GeV/n 
for both species. This is similar to our results for protons but higher (400 GV) than our results 
for helium. Results from ATIC-2 (18) implied that protons and helium nuclei have different 
energy spectra, although the results suffered from unclear systematic uncertainties and there 
were differences with respect to previously reported ATIC- 1 (33) data. 
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